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Abstract. We consider Casimir force acting on a three dimensional rectan- 
gular piston due to a massive scalar field subject to periodic, Dirichlet and 
Neumann boundary conditions. Exponential cut-off method is used to derive 
the Casimir energy in the interior region and the exterior region separated by 
the piston. It is shown that the divergent term of the Casimir force acting 
on the piston due to the interior region cancels with that due to the exterior 
region, thus render a finite well-defined Casimir force acting on the piston. 
Explicit expressions for the total Casimir force acting on the piston is derived, 
which show that the Casimir force is always attractive for all the different 
boundary conditions considered. As a function of a — the distance from the 
piston to the opposite wall, it is found that the magnitude of the Casimir 
force behaves like I/a 4 when a — > 0+ and decays exponentially when a — > oo. 
Moreover, the magnitude of the Casimir force is always a decreasing function 
of a. On the other hand, passing from massless to massive, we find that the 
effect of the mass is insignificant when a is small, but the magnitude of the 
force is decreased for large a in the massive case. 

1. Introduction 

Casimir effect associated with piston geometry has attracted considerable inter- 
est since its introduction by Cavalcanti pQ a few years ago. The main attraction 
of Casimir piston is that such a geometric setup can resolve a serious divergent 
problem that plagues the Casimir calculations. In the conventional calculations 
of Casimir force inside a confined region such as a rectangular box, the nontrivial 
contribution of vacuum energy outside the box has been ignored, and the surface 
divergent terms which depends on the dimensions and geometry are discarded to 
obtain finite results under the pretext of some regularization schemes. It has been 
shown that such surface divergence cannot be removed by renormalization of phys- 
ical parameters of the theory [2j [3]. Calvacanti [I] considered a new geometric 
configuration - a two dimensional rectangular piston, and showed that for a mass- 
less scalar field subject to Dirichlet boundary conditions, the divergent part of the 
Casimir force acting on the piston due to the interior and exterior regions cancel 
with each other and the net result is a finite attractive Casimir force without the 
surface divergence. This result has stimulated an interest in studying the Casimir 
force on pistons with different geometric setups and boundary conditions. In par- 
ticular, Hertzberg et al [3] have studied the Casimir effect for electromagnetic fields 
in three dimensional rectangular pistons with perfectly conducting walls. It was 
found that the Casimir force is always attractive. This work was generalized in 
[5] where finite temperature effect was taken into account and pistons with general 
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cross sections were considered. Around the same time, Marachevsky also carried 
out a similar investigation using a different approach El [8] . In [9] , Edery studied 
the case of a massless scalar field with Dirichlet boundary conditions for three di- 
mensional rectangular pistons. He also found that the force on the piston is always 
attractive, in contrast to the fact that the regularized Casimir force for a three 
dimensional rectangular cavity for massless scalar field with Dirichlet boundary 
conditions can be attractive or repulsive depending on the relative size of the cav- 
ity. This work was later generalized in [TU] to massless scalar fields with Dirichlet 
and Neumann boundary conditions in any dimensions. In all the above scenar- 
ios, the Casimir force was found to be attractive. In [11] , Barton found out that 
when weakly reflecting dielectric materials are used, the Casimir force on the three 
dimensional piston can become repulsive when the plate separation is sufficiently 
large. On the other hand, the work of Zhai and Li [12] showed that in the case of 
mixed boundary conditions (one plate with Dirichlet boundary conditions and one 
with Neumann boundary conditions), the Casimir force on a rectangular piston in 
one, two and three dimensions is always repulsive. More recently, Casimir effect 
has been investigated for electromagnetic fields with perfect magnetic conditions in 
rectangular pistons of arbitrary dimensions [13j and for massless scalar fields with 
Dirichlet boundary conditions in pistons inside space-time with extra compactificd 
dimensions [T3] . 

To the best of our knowledge, no work has been carried out on the Casimir effect 
for massive scalar fields in the piston setting. In this work, we consider a massive 
scalar field in a three dimensional rectangular piston subject to periodic boundary 
conditions, Dirichlet boundary conditions as well as Neumann boundary conditions. 
We use exponential cut-off method to compute the cut-off dependent Casimir energy 
inside and outside the piston. We show that the sum of the Casimir energies has a 
divergent part that does not depend on the position of the piston, and therefore, the 
force acting on the piston is finite, without any regularization. Explicit formulas for 
the Casimir force are derived, which can be written as an infinite convergent sums 
of Bessel functions. From the formulas, it is easy to deduce that the Casimir force 
is always attractive under all the different boundary conditions. We also proved a 
stronger result: as a function of a — the separation between the piston and the 
opposite wall, the magnitude of the Casimir force behaves like 1/a 4 when a — * + 
and it decays exponentially when a — + 0. Moreover, the magnitude of the force 
is always decreasing from infinity to 0. The mass effect to the Casimir force is 
considered and it is found that the mass effect is significant only when a is large. 
Some numerical simulations have been carried out. 

2. Cancelation of divergence in three dimensional Casimir piston 

Consider a 3-dimensional rectangular piston in the form of a rectangular cavity 
[0, Li] x [0, L 2 ] x [0, L 3 ] separated by a plane x\ = a (the piston) into two regions: the 
interior region [0, a] x [0, L 2 ] x [0, L 3 ] and the exterior region [a, L\] x [0, L 2 ] x [0, L 3 ] 
(see Figure 1). At the end we are going to let L\ — > so that the exterior region 
becomes an open region. We want to compute the Casimir force acting on the 
piston due to a massive (to > 0) scalar field <fi(x., t) satisfying the Klein-Gordon 
equation in Minkowski space-time: 



(2.1) 





THREE DIMENSIONAL CASIMIR PISTON FOR MASSIVE SCALAR FIELDS 



3 





Figure 1: The three dimensional rectangular pistons 



where 



□ 
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dt 2 dx 2 ' 



The periodic boundary condition (be), Dirichlet be and Neumann be will be con- 
sidered. 

The Casimir energy of the piston system is the sum of the Casimir energy of 
the internal region _E cas (et, L 2 , L3; m) and the Casimir energy of the external region 
E cas (Li — a, L 2 , L 3 ; m). The Casimir force acting on the piston is then obtained by 

(2.2) F Cas = lim --^( E Ca s(a,L 2 ,L 3 ;m) + E Cas (Li ~ a, L 2l L 3 ;m)\ 

In the periodic be case, the eigenmodes of the (d + 1) - dimensional field c/)(x, t) 
in a d-dimensional rectangular cavity [0, L{\ x . . . x [0, Ld] satisfying (12. 1| (with 3 
replaced by d) are given by 



2tt 



+ m l 



k £ Z d : 



and the Casimir energy is defined as the divergent sum 
(2.3) E^ s (L 1: ...,L d ;m) ' 



2 E 

kez d 



For the Dirichlet be and Neumann be cases, the eigenmodes are respectively 

+ m 2 , keN d ; 
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Therefore when d = 3 [15] , 
(2.4) 



E^ s N (L 1} L 2 ,L 3 ;m) =-<j ££ as (2L 1( 2L 2 , 2L 3 ; m) 



1 



T Pcb8(2£i, 2£ 2 ; m) + Eg^pLt, 2L 3 ; m) + E^(2L 2 , 2L 3 ; m) 
+ ££ as (2£i; m) + E^ s (2L 2 ;m) + E^ S (2L 3 ; m) T ~m I. 



Different regularization schemes have been employed to define a regularized Casimir 
energy from (|2.3[) . We adopt here the exponential cut-off method which allows us 
to retain the divergent terms. Define the A-dependent Casimir energy by 



(2.5) 



,P„-Au 



In the Appendix A, we show that 
(2.6) 



-EtLs (A;£i, ...,L d ;m) =A£ as , div (A;£i, . . .,L d ;m) + E^ s reg (ii, . . .,L d ;m) + 0(A), 

where for rf = 1,2 and 3, the A — > + divergent term S^as div > • • • > ^«*; TO ) 

r r , 

Cas, regv 



and the regularized Casimir energy -E^as rcE (^i 5 • ■ ■ , L d ; m) are given respectively 



by 



E, 



Cas, div 



(A;Li;m ) =— A — logA; 

27T 47T 



(2.7) 
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(2.9) 



r r t L\L 2 L 3 4 / 1 \ L\L 2 L 3 2 

- ~^T m ("2 +7 + logm-log2j -^-^m 2 



^Cas, rcgC^lj ^2, ^3! TO,) 



x S ( S [ fciLj ] 2 J ^ 2 ( m A 

kez3\{o} \t=i / 



\ i=l 



Here 7 is the Eulcr constant and K v {z) is the modified Bessel function. From 
these formulas and (|2.4p . we can compute the cut-off dependent Casimir energies 
for massive scalar fields inside three dimensional rectangular cavities under various 
boundary conditions. Conventionally, a finite Casimir energy is defined by taking 
A — > + limit in the regular term of the Casimir energy. We would like to comment 
that the result will agree with the result derived by zeta regularization method 
(see Appendix B) if one defines the normalization constant fi appearing in the zeta 
regularization method as e -7 . 

Now we analyze the divergent term of the Casimir energy. In the periodic be case, 
we observe that the divergent term depends linearly on the volume V — L\L 2 L 3 
of the rectangular cavity. However, for the cases of Dirichlet be and Neumann be, 
(|2.4p shows that the divergent term can be written in the form 

(2.10) 

E Ca^div( X '> L ^ L ^ L 3'> m ) 

=£>y(A; m)LiL 2 L 3 ± £> S (A; m)(L x L 2 + L 2 L 3 + LxL z ) + D L (A; m)(L 1 + L 2 + L 3 ). 



In other words, besides the bulk divergence Dv(^',m)L 1 L 2 L 3 which only depends 
on the volume, there are also surface divergence D5(A;m)(X 1 L 2 + L 2 L 3 + LxL 3 ) 
and divergence due to corners £>l(A; m)(Li + L 2 + L 3 ). 

For the piston geometry, it is easy to see from (|2 . 1 0|) that the divergent term for 
the sum of the Casimir energy of the interior region and the exterior region given by 
Eq(^ (A; a, L 2 , L 3 ; m) + E^[^ N (X; Li — a, L 2 , L 3 ;m) depends only on Li,L 2 ,L 3 , 
but not on the position of the piston xi = a. In other words, the divergent parts of 
the Casimir energies in the interior and exterior regions contribute Casimir force of 
same magnitude but opposite signs to the piston, and thus cancel with each other. 
This renders a finite quantity for the Casimir force acting on the piston given by 
(|2.2p . which can be rewritten as 



(2.11) 

F "- /N - ^|(C 8 (^ = 0; a, L 2 , L 3 ; m) + 2^£(A = 0; L x - a, L 2 , L 3 



3. Analysis of the Casimir force 

From (|2.7|) , (|2.8[) and (|2.9p . one observes that the regularized Casimir energy 
i^Cas, reg (A = 0; L\, . . . , Ld\ m), d = 1, 2, 3 inside a rectangular cavity can be written 
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d+1 



as a term 23^(^2, ■ • ■ , Ld] m)L\ which depends linearly on L\, plus a Bessel series 

2 feiez\{o} \i=i j 



x Kd+i m» 



d 

[kiLi 



Notice that the summation over k e Z d \ {0} in (|2~7)) . ([278]> and (|2~9|l is de- 
composed into summation over fci = 0, (fej, . . . , k d ) £ Z d_1 \ {0} and summa- 
tion over k\ ^ 0, (k 2 , • • • , kd) G Z d_1 . The term that involves summation over 
A'i = 0,(k 2 , ■ ■ ■ ,k d ) G Z d_1 \ {0} depends on L\ linearly and is therefore com- 
bined into the term 53^(^2, ■ • • , Ld] m)L\. As in the case of the divergent terms, 
for periodic be as well as Dirichlet be and Neumann be, the sum of the terms 
Q3d(£ 2 , ■ ■ • , Ld] m)L\ of the interior region and exterior region does not depend on 
the position of the piston and thus the corresponding Casimir force cancels with 
each other. In other words, the Casimir force acting on the three dimensional piston 
is given by 

(3.1) 

F c as ( a \ L 2, L 3 ;m) = - lim — \ R 3 (a, L 2l L 3 ; m) + R 3 (Li - a, L 2l L 3 ; m) \ 

Li^oa da V ) 



in the case of periodic be; whereas for Dirichlet and Neumann be, (|2.4|) shows that 

(3.2) F^ s N (a;L 2 ,L 3 ;m) 

= -l lim ^-\R 3 (2a,2L 2 ;2L 3 ;m) + R 3 (2L 1 -2a,2L 2 ,2L 3 ;m) 
o Li^oo oa l 

T R 2 (2a, 2L 2 ; m) T # 2 (2Li - 2a; 2L 2 ; m) T R 2 {2a, 2L 3 ; m) 
T i? 2 (2£i - 2a,2L 3 ;m) + R 1 {2a;m) + i?i(2L x - 2a;m)|. 
Using the formula ([16], #3.478, no. 4), 

(3.3) / t ±v - x exp (-at - dt = 2 ( -] * K v {2yQ3), 



a 



we have 
(3.4) 

R d{Ll ,...,L d] m) = -%^ [°°t-*P Yl expf-^-l^M 

(fc 2 ,...,fc d )ez d - 1 

From here, we see that R d (Li, . . . , Ld] m) decays exponentially as L\ — ► oo. Taking 
derivative with respec to a and letting L\ — > oo, we have 

d 

lim -—R d (L 1 -a,...,L d ;m)=0. 

L1-+00 oa 

In other words, there is no contribution to the Casimir force on the piston from the 
terms Rd{L\ — a, L 2 , . . . , Ld] m) of the exterior region. Therefore, the net contri- 
bution to the Casimir force on the piston are due to the terms Rd(a, L 2l . . . , Ld', m) 
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of the interior region. More precisely, 



(3.5) 
and 



d 



F Gas( a j L 2, L a ;m) = —Q-Rzia, L 2 , L 3 ;m), 



lid c) 

(3-6) ^caf( a ; m) =-< —^Ra(2a, 2L 2 , 2L 3 ; m) ± ^i?a(2a; 2L 2 ; m) 

± t^M^ 2L 3 ; m) - — i?x (2a; m) I . 
Now applying the Jacobi inversion formula 

oo i — oo 
k— — oo /c— — oo 

to the summation over (fc 2 , . . . , kd) £ in (|3.4p . we find that 

(3.8) 



Rd(a 1 L 2 , ...,L d ;m) 



8rr Jo 



poo / d 

I r 2 ex P ~ tm2 ~ 4i7r2 E 



feiSZ\{0} 

{k 2 ,...,k d )& 



hi 
Li 



At 



[aki] 2 dt 



1 °° 1 



r 27rfc," 1 " 



E 



^ fci\ 

fel=l (fc 2 ,...,fc d )eZ' 1 -i ^ i=2 

Eq. p.5p then shows that in the periodic be case, 



K\ ak\ 



\ 



E 

z=2 



u 



1 °° 1 

"tt E E Y x 

Aii=l (fc 2 ,fc 3 )GZ 2 



\ 



E 

i=2 



27rfc, 



I ah 



\ 



E 

i=2 



27rfc, 



On the other hand, using the fact that if f(k 2 ,k 3 ) is a totally even function, i.e. 
f(±k 2 ,±k 3 ) = f(k 2 ,k 3 ), then 



3 £ t E t E 

I,(fc2,fc3)ez 2 fc 2 ez,fc 3 =o fc 2 =o,fc 3 ez fc 2 =o,fc 3 =o 
Eqs. (|3.6[) and (|3.8p give the expression 



J2 \f(k 2 ,h)= f( k * 

(fe 2 ,fc 3 )eN 2 /(Nu{o}) 2 
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for the Casimir force in the Dirichlet be and Neumann be cases. Now using the 
formula [[16 , #8.486, no. 12] 

(3.9) Tz K ^ z) = ~~ z Kv{ ^ ~ K v-i( z )> 

we obtain finally the explicit expression for the Casimir force in the periodic be, 

Dirichlet be and Neumann be cases: 

(3.10) 



^ca S ( a ; L 2,£3;m) 
1 00 1 

— V V — 

ira ^ ^ k\\ 

k 1= i(k 2 .k 3 )ei 2 \ 

-It E f- 



E 



K\ aki 



\ 



E 



2itki 



3 r 



ki=i (fe 2 ,fe 3 )ez 2 



E 

i=2 



2nki 



1 2 N 



Li 



K ah 



\ 



3 r 



E 

i=2 



2nki 



~] 2 



(3.11) 



F caf r (a;L 2 ,L 3 ;m) 



1 oo 

— E 



1 



2lTU ^— E \ 

fci = l (fc 2 ,fc 3 )GN 2 /(NU{0}) 2 \ 



+E 



i=2 



17 



Kx 2ak x 



i 



+ E 



(=2 



7Tfc,; 



E E 



3 r 



fei = l (fc 2l fc 3 )eN 2 /(Nu{0}) 2 



E 

i=2 



\ 



E 

i=2 











[Li. 





These formulas are useful for studying the large-m and large-a behavior of the 
Casimir force. More precisely, since the modified Bessel function K v {z) decays 
exponentially as z — > oo, with leading term 



K„{z) 



Eqs. (|3.10[) and (|3 . 1 1 [) show that for fixed m, L 2 and L$, in the large a (a 3> 1) 

Cas 



limit, the Casimir force FniF {a; L 2l L&m) decays exponentially with leading 



terms being 
Fc m (a;L 2 ,L 3 ;m) 



-m\j — e , 
■na 

' 2V na 



1 / 1 

2 V 



2 2 \ 3/4 

7T TT ' 

L 2 L 3 



exp —2a* m? + 



respectively. In other words, for fixed L 2 and L3, when the piston is very far 
away from the opposite wall, the Casimir force is an exponentially weak attrac- 
tive force. In the massless m — > + limit, the Casimir force F^(a; L 2 , L3; 0) for 
Dirichlet be still decays exponentially with respect to a. However, the Casimir force 




^Cas( a ! -^2, L3; 0) and F^ s (a; L2, L3; 0) for periodic be and Neumann be becomes 
polynomially decay with order 1/a 2 . 

Another advantage of the explicit formulas (|3.10|) and ()3. 1 ljl is that since the 
Bessel function K u (z) is positive for all z > 0, Eqs. (|3.10p and (|3 . 1 1[) show im- 
mediately that the sign of the Casimir force F^ a f^ N (a; L2, L3; m) for all the three 
boundary conditions is negative, and therefore they are all attractive forces. This 
extends the known results about the massless scalar fields [TJ [9l [10] . Extension of 
this result to massive scalar fields in arbitrary dimensions will be considered in a 
future work. Here we would like to remark that as is shown in Appendix B, the 
Casimir force acting on the piston due to the regularized Casimir energy of the 
interior region alone can be either attractive or repulsive for any of the boundary 
conditions. It is the total effect of the Casimir force from both the interior and 
exterior regions that gives a net finite force which is attractive. 

For the behavior of the Casimir force when the separation distance a is small 
compared to L2 and L3 , we apply again the formula (|3.7[) to the summation over k\ 
term of (|3.4[) . The detailed derivation is given in the Appendix C. Here we present 
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the end results: 
(3.12) 



L 2 L 3 2 



27r 2 a 2 



00 1 L 2 L 00 1 

E-T T K 2 (afnk 1 ) — - „ m 3 > — 
A:? v ; 2n 2 a ^ k x 

<ci = l 1 fci = l 1 



(amk\ ) 



L 2 L. 
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2V2ttL 2 L 3 



(fc 2 ,fc 3 )6Z 2 \{0} \i=2 / 



E E *? 1 

fe 1= i (fc 2 ,fc 3 )ez 2 \{o} VX^I^*] 2 




2 t 2 \ 1/4 



for the periodic be case; and 
(3.13) 



L 2 L 3 2 



OO ^ 

^ -rjK 2 (2amki) 

fe 1= i 1 

(£2 + £3) to 3/2 ^ 3/2 
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Air 2 a 



oc _^ 

V —K^amh) 



fei = l 



8^ 3 / 2 
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E k 1 3/2 K 3/2 (2amk 1 )± 



fci=i 



(£ 2 + £3) m 5 / 2 ^ t _ 1/2 
4tt 3 / 2 



1/2" E fc i ^5/2(2amfci) 



fe,=i 



1 m v ■> 1 m x -> 
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Rtt ^ ^ — ^ Zfi 47T — ' 
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^3 fe 3 =i 
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16tt 2 



E^] 2 K2 \^ 
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E E *? 

fci=l (fc 2 ,fe 3 )6Z 2 \{0} 



#1/3 



ill 



TT 2 k 2 



for the Dirichlet be and Neumann be cases. The last term on the right hand side 
of (|3.12p and the last three terms on the right hand side of (|3.13p vanish to zero 
exponentially fast when a — > + . The third term on the right hand side of (|3.12p and 
the seventh, eighth and ninth terms on the right hand side of (13.13[) are independent 
of a. For the remaining terms, their behaviors as a — > + are not obvious. Using 
the power series expansion of the Bessel function K v (z) naively do not give us the 
correct asymptotic behavior due to the summation over h\. In Appendix C, we 
derive the correct asymptotic behavior for these terms when a — > + . The results 
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are: for fixed to, L 2l L 3 , as a — > + , 

, T T v tt 2 L 2 L 3 L 2 L 3 m 2 7r 2 m 4 , n 

^Cas(a; £2, £3; m) -r-j- + — 2 + ——L 2 L 3 \oga + 0(a u ) 

3Uaf 24a z 32 

pD/jv, T j- x 7r 2 L 2 L 3 C(3)(^2+-^3) 7T . L 2 L 3 m 2 m 2 (L 2 +L 3 ) 

^Cas W -^2, ^3! TO.; — — 4 ± — 5 — — H 2 T 

480a 4 167rcr 96er 9oer 167ra 

2 4 2 
7T TO T r , 7TTO , n , 

+ -^—^2-^3 log a — loga + 0(a u ). 

32 16 

In particular, as a — > + , the Casimir force becomes very negative (attractive) with 
leading term of order 1/a 4 . Notice that this leading term is independent of the 
mass to. We also observe that the first three leading terms of the Casimir force 
F® as (a; L 2 ,L 3 ; to) for Dirichlet be has exactly the same form as that of the massless 
case derived in [S]. In fact, by taking the limit to — > + for F^ as (a; L 2 , L 3l m) in 
(|3.13p , we obtain the same expression as that derived in [5] when a is small. 

From the above discussion, we conclude that the effect of mass is insignificant 
when a is small, but with nonzero mass, the Casimir force will be drastically re- 
duced. 

We would also like to mention that the expressions (|3.12[) and (|3.13[) can also 
be used to analyze the behavior of the Casimir force when L 2 and L 3 are large 
(parallel plate geometry). In this case, we find that the Casimir pressure (force per 
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0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 



a a 

Figure 4: These figures show the dependence of the Casimir force FQ as (a; L2, 1/3; m) (pe- 
riodic be) on a when L2 and L3 are varied. On the left hand side, L3 is fixed to equal to 
1 and 1/2 = 1, 2, 5, 10. On the right hand side, L2 is set to equal to L3. 



L = 1 ; m = 1 m = 1 

-10-*, , , , , , I 1 1 




° 02 M 0.6 0.8 1 0.2 0.4 0.6 0.8 1 



Figure 5: These figures show the dependence of the Casimir force F Cas (a; L2, L3; m) 
(Dirichlet be) on a when L2 and L3 are varied. On the left hand side, L3 is fixed to 
equal to 1 and L2 = 1, 2, 5, 10. On the right hand side, L2 is set to equal to L3. 



unit area) acting on one of the plates is 

D p , . ,. F Cas (a;L 2 ,L 3 ;m) 

P Ce.sW m ) : = r p — 



2n 2 a 2 ^ fc? y ' 2ir 2 a ^ k x 
fc 1= i 1 fci=i 

DIN , m 1 m 1 

p cL («;™) =^2T2 2^ ^K 2 (^mki) - ^ 2^ J-K^amh). 
fci=i 1 fci=i 



In the massless limit m — > + , these become Pc as (a; 0) = — 7r 2 / (30a 4 ) and P Cas (a; 0) = 
— 7r 2 /(480a 4 ) respectively. 
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L = 1; m = 1 
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Figure 6: These figures show the dependence of the Casimir force Fq^w, L2, L3; m) (Neu- 
mann be) on a when L2 and L3 are varied. On the left hand side, L3 is fixed to equal to 
1 and L2 = 1, 2, 5, 10. On the right hand side, L2 is set to equal to L3. 



In Figure 2, 3, 4, 5, 6, 7 we show graphically the behavior of the Casimir force 
Fq[^ N (a; L 2 , L 3 ; m). From (|3.10[) and (|3.11[) . it is easy to verify that 

FZ D/N (Aa; XL 2 , AL 3 ; 2) = X^F^ia; L 2 , L 3 ; m). 

Therefore, we have the freedom of fixing one of the variables which does not affect 
the shape (up to scaling) of the graphs. When one is not concerned with the 
dependence of the Casimir force on mass, one can let m = 1. From Figs. 2, 4, 5, 6, 
we see that fixing m, L 2 and L3, the Casimir force always increase (as a function 
of a) from — 00 (when a — > + ) to (when a — > —00). In fact, by differentiating 
(I3.10|) and p. lip again with respect to a and using Q3.9p . we can conclude that 
the Casimir force Fq[J^ (a; L 2 , L^.m) is always an increasing function of a, and 
therefore the small a and large a behaviors show that it always increases from —00 
(when a — ► + ) to (when a — ► —00). 

Fig. 3 shows the dependence of the Casimir force on mass m. We see that the 
mass effect is more significant when a is large. On the other hand, we will tend to 
make conclusion that the Casimir force is an increasing function of m from Fig. 3. 
However, Fig. 7 invalidates this conjecture. 

In Figure 8, we demonstrate the Casimir force on the piston, the regularized 
Casimir force and reduced regularized Casimir force (see Appendix B) due to the 
interior and exterior regions on a single graph, for the periodic be and Dirichlct 
be cases. The graphs show that the Casimir force due to the interior region and 
exterior exterior region can be cither attractive or repulsive. But the net effect on 
the Casimir piston is attractive. 

In the above discussion, we are considering the case where the limit L\ — > 00 
has been taken. Now we briefly discuss the case where the limit L± — > 00 is not 
taken, i.e. the piston is assumed to be placed in a fixed closed rectangular box. To 
distinguish with the cases above, we denote the resulting Casimir force on the piston 
as FqI^^ N (a; L 2 , L 3 ;m). From the discussion in the beginning of this section, 
Fq^ 3 (a; L 2 , Lz^m) is given by the right hand side of (|3.1|) and (|3.2p without 
the limit limz, 1 ->. 00 - Therefore, in terms of the Casimir force F^[^ N (a; L 2 , L^;m) 
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a = 1; L = 1; m = 1; periodic be a = 1; L = 1; m = 1, Dirichlet be 




Figure 8: These figures show that despite the fact that the Casimir force acting on the 
piston is always attractive, the regularized Casimir force due to the interior region and 
exterior regions can be either attractive or repulsive. 



computed for the case L\ — > oo, we have 

Fcaf /N ( a > L 2> L 3, m ) = FcU? (a; L 2 ,L 3 ; to) - F^' N ' {L x - a; L 2l L 3 ; to) 

Obviously F^[^ N (a; L 2 , £3; to) = when a = Li/2, i.e. the piston is in an 
equilibrium position when it is placed in the middle of the box. However, if it is 
placed closer to the left end, i.e. a < L± — a, then since Fq^P (a; L 2 , L^m) is 
more negative than F^[f^ N (Li — a; L 2 , L3; to), we find a net attractive force pulling 
the piston to the left end. In other words, the Casimir force always tends to move 
the piston to the nearest wall. In nanotcchnology, this is known as the stiction 
effect and may have undesirable consequences to the functionality of nano-devices. 

4. Conclusion 

We have carried out the study of the Casimir effect for massive scalar fields in 
three dimensions subject to different boundary conditions in the piston setting. 
The methods we used here can be easily generalized to show that for pistons with 
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arbitrary cross sections in any dimensions, the Casimir force acting on the piston 
is always divergence free. For the cases considered here, the Casimir force is found 
to be always attractive. It would be interesting to investigate the possible repulsive 
Casimir force due to massive fields in the piston setting. Another thing we have 
not addressed here is the thermal effect, which we plan to carry out in the future. 
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Appendix A. Cut-off dependent Casimir energy for massive scalar 

FIELDS IN RECTANGULAR CAVITIES 

The cut-off dependent Casimir energy E^ as (X; L\, . , . , L d ; to) (|2.5p for a massive 
scalar field with periodic be inside a <i-dimensional rectangular cavity can be written 
as 

1 d 

(A.l) E^ as (\;Li, ...,L d ;m) = --—K(X;Li, . . .,L d ;m), 



where 



K(X;L 1 , ...,L d ;m) = 



For a massless (to = 0) scalar field, the sum (|A.1[) has been computed in [T7] using 
Euler-Maclaurin summation formula. Here we want to compute this kind of sum 
by a different method, making full advantage of the analytic structure of Epstein 
zeta function. Using inverse Mellin transform of exponential function 



1 

2ni 



C0+2OO 

T{w)z- W dw, c > 0, 

cq — ioc 



we find that 
(A.2) 

K{\;L 1 ,...,L d ) = — \- w T(w)Z d (—;—,..., —;m ) dw, c > d, 

2m J Co _ ioo V 2 L i L d J 

where Z d (s; oi, . . . , a d \ to) is the inhomogeneous Epstein zeta function defined by 
the series 



Z d (s; ai,...,a d ;m)= [^h] 2 + 

kez^ \i=i 
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when Re s > d/2. Zd{s; ai, . . . , ad] to) has a meromorphic continuation to C given 

by (see e.g. USES EDI E5) 

(A.3) 



Z d (s;ai, ...,a d ;m) = 



2tt s 



r( s ) 



E E 

keZ d \{0} \i=l 



if. 



where K v (z) is the modified Bessel function. The second term is an analytic func- 
tion of s. The first term shows that when d is odd, Zd(s] a\, . . . , ad] m) has simple 
poles at s = i — j,j = 0, 1, 2, . . . with residues 

(— IV 7r^ to 2j 
(A.4) HfiB a= ^_ J .Z < ,(a;o 1 ,... J o d ;m) = — — — j ± _ . 

J ' lU=i a « V2 

When rf is even, Zd(s; a\, ■ ■ ■ , ad',m) only has poles at s = | — j where j = 
0, 1, . . . , with the same residue (IA.4|1 . By moving the contour of integration 
from the line Re w = Cq, c > d to the line Re w = —2 + s, we obtain for d — 1, 2, 3, 



K{\;L x ;m) =— A -1 + 



Xm 2 Li Am Ij /. , „ 5 

— log A H — I log m — log 2 — ^ + 7 



2Am 



^fc- 1 ifi(mfcL 1 ) + 0(A 2 ); 



fe=i 



K(X]L 1 ,L 2 ]m) 



Am 3 + 



L1L2 2 L1L2 2 , -^1-^5 

: — — A — m H — - — 

27T 47T D7T 

x E [E^n ^3/2 1"'* 

keZ 2 \{0} \i=l / 



TIP 



Am 3 / 2 



\ E t^-^*] 



0(A 2 ); 



and 



K(X;Li,L2,L 3 ;m) 



L1L2L3 



A" 



L1L2L3 2\— 1 L1L2L3 

-to A 



4tt 2 



16tt 2 



L1L2L3 4 , 



16tt 2 



to A h 7 + log m — log 2 



to A log A 
L1L2L'. 



2tt 2 



Ato 2 



\ E 



0(A 2 ), 



x E^M A '-> I 

k6Z 3 \{0} \i=l / 

where 7 is the Euler constant. Using (|A.1[) . we find that the cut-off dependent 
Casimir energy E^ as (A; L\, . . . ; L 3 ; to) for d = 1, 2 and 3 are given respectively by 



^(A;£i;m) 



2tt 



lA- 2 



m 2 Li m 2 Li 
■log A 



4tt 



4tt 



log to — log 2 — - + 7 



- Vfc- 1 X 1 (TOfcL 1 )+0(A); 

7T ^— ' 



fc=l 



THREE DIMENSIONAL CASIMIR PISTON FOR MASSIVE SCALAR FIELDS 



17 



P p <\-T T -ml - LlZ/2 \-3 LlL2 ^3 L l £ 2 3/2 

^Cas( A ) -^l) -^2, m; — — — A — — m i==m 



2tt 



12tt 

' 2 



'87T 3 



x Etw 

keZ 2 \{0} \i=l 



#3/2 I m, 



, ^ [kiLi 



0(A); 



and 



-^Cas (A;Li,L 2 , J L 3 ;m) 



3L1L2L3 4 L1L2L3 2\-2 L1L2L3 



2vr 2 



-A 



L1L2L 



3_4 



32tt 2 



8tt 2 
1 



m^A~^ + 



32rr 2 



-m log A 



m — — + 7 + log m — log 2 — 



L1L2L3 
4tt 2 ' 



keZ3\{0} \i=l / 



^[fc 4 L,] 2 ] + 0(A). 



One can show that in the massless limit (m = 0), these results coincide with the 
corresponding results obtained in [17l [TO] . 



Appendix B. Casimir effect for massive scalar fields in rectangular 

CAVITIES BY ZETA REGULARIZATION METHOD, AND THE 
ATTRACTIVE/REPULSIVE NATURE OF THE CASIMIR FORCE 

Using zeta regularization method, the Casimir energy Eq (L\, . . . , L^; m) (|2.3|) 
for a massive scalar field with periodic be inside a o?-dimensional rectangular cavity 
is defined as [25] 



£&.(£i, ■ • • , L d ; m) := 1 lim Up f + £ j ^ 



where Cp( s ) is the zeta function 



,P\-2s 



2tt 



2tt 



= ^ s Z d [s;— ;m 



and is a normalization constant. Using (|A.3[) . we find that E^^Lx, . . . ,Ld\m) 
depends on the normalization constant /1 if and only is d is odd. For d = 1,2,3, 
-^Cas( L i> ■ ■ ■ > L d\ m) is given explicitly by 



L L I 1 \ 

E c as ( L i; m ) = ^m 2 logn- ( logm-log2- - J - — Y k ~ lK i( mkL i)> 



Ec as {Li,L 2 ]m) = - 



L1L2 3 il-C/2 3/2 

12tt V8t73 

/ 2 \ " 



x J2 Efe £ «n ^/2 

keZ 2 \{0} \i=l / 



\ E 

\ »=i 
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L1L2-L3 4 L1L2L3 4 



32tt 2 



L1L2L 



-m log/i 



32tt 2 



logm - log2 - - 



3^2 

m 



4rr- "' ^ . 

k£Z3\{0} \i=l 



Ki 



\ E t^*^ 



Compare to the result obtained by cut-off regularization method, we find that 
the zeta regularized Casimir energy Eq^Li, . . . , La; to) will agree with the A- 
independent part of the cut-off dependent Casimir energy Eq^X; L±, . . . ,Ld',m) 
if and only if we take fi = e -7 , where 7 is the Euler constant. Using this prescrip- 
tion, we find that the zeta regularized Casimir force acting in the X\ direction for 
d = 3 is 

(B.l) F 1 % s (L l7 L 2 , L 3 ; to) = -A-^Lt, L 2 , L 3 ;m) 



L 2 L 



2-^3 4 
to 



32tt 2 



logm - log 2 - - + 7 



2-^3 2 
TO 



4tt 2 



2 TO, 



(/c 2 ,/c 3 )£Z2\{0} \i=2 



\ E i^^' 

\ i=2 



/ 



L1L2L3 

4tt 2 ' 



2 E^ L d 2 ^ K 



V 



fciSZ\{0} \<=1 

(fe2,fe 3 )ez 2 



\ ,-=i 



\ 

I 

/ 



As is shown in Section 3, the last term in (jB.ip is the net contribution to the 
Casimir force acting on the Casimir piston. For any fixed m,L 2 >^ 3 j h is always 
negative and increase from — 00 (when L\ — > + ) to (when Lj — > 00). On the 
other hand, the first and second term in (|B. 1[) are independent of L±. Moreover, 



the second term is always positive and by taking m sufficiently small, the first term 
is also positive. Therefore, for fixed L 2l L 3l when m is sufficiently small and L\ is 
sufficiently large, the zeta regularized Casimir force Ff Cas (Li, L2, L3; to) becomes 
repulsive. Similar reasonings can be used to show that the zeta regularized Casimir 
force F^^ s (Li, L 2 , £3; to) acting in the x% direction of a rectangular cavity due 
to a massive scalar field with Dirichlet and Neumann boundary conditions can be 
either attractive or repulsive depending on the values of L\, L 2 , L3 and m. 

We would also like to remark that in some of the regularization schemes, the 
Casimir force on a wall in a rectangular cavity is defined by subtracting the Casimir 
force in the absence of boundary. In this case, the Casimir force is equal to the sum 
of the last two terms in (jB.ip . i.e. the first term is absent. We call the Casimir 
force defined in this way the reduced Casimir force. Notice that in the massless 
case, there is no difference between the Casimir force and reduced Casimir force in 
a rectangular cavity. On the other hand, it is easy to see that the argument given 
above shows that the reduced Casimir force in a rectangular cavity can be positive 
or negative, under different boundary conditions. In Figure 8, we demonstrate the 
Casimir force on the piston, the Casimir force (and reduced Casimir force) due to 
the interior and exterior regions on a single graph, for the periodic be and Dirichlet 
be cases. 
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Appendix C. An alternative expression for the Casimir force when a 

IS SMALL 

To study the behavior of the Casimir force when a is small, we first split the 
summation over (fc 2 , . . . , kd) G Z d_1 in (|3.4|) into the term (fc 2) . . . , kd) = and the 
summation over (fc 2 , . . . , kd) G Z d_1 \ {0}. Then applying the identity (|3.7|) to the 
summation over fci G Z \ {0}, we find that 

(C.l) 

Rd(a, L 2 , ■ ■ ■ , L d ; m) 

z ^ " (fe2,-,fe<j)ez <i - 1 \{0} \ i=2 / 

- /„ «- ^ e - (-< r + -vj - 5 § j *■ 

(fc 2 ,...,fc d )ez d - 1 \{o} 
Taking derivative with respect to a gives 
(C.2) 

d 

- -^-Rd(a,L2, ...,L d ;m) 
oa 

= lLl - 2 ' / t'— \ exp -tm 2 - dt 

2*+^ Jo ^ v 4t ; 

z ^ (fc 2 ,...,fc ti )ez<i- 1 \{0} \ i=2 / 

+ 2 it" 3 2 / *-* E fe2 

(fe 2 ,...,fc d )ez £i - 1 \{o} 

The expressions for the Casimir force i*oae (a; Li, L 2 , L3; to) can then be ob- 
tained by using (|3.3[) . (|3.5p and (|3.6p . The results are given by (|3.12[) and (|3.13|) . 

The last term on the right hand side of (|C2[) tends to zero exponentially fast 
when a — > + . The third term is independent of a. For the first two terms, their 
asymptotic behaviors when a — ► + are not obvious. Let us call the sum of these 
two terms as T(a;L 2 , . . . ,L d ;m). Here we derive its asymptotic behavior when 
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a — > + . From the derivation above, we have 

cr, t ^ d i a Ui-^ L i f°° - d + 3 ^ 

S(o; L 2 ,..., L d ; m) =— I d+1 / t * \ exp 



j fc 2 a 2 
4f 



nt 2 ^ d m f d+i 



where *P(s;a;m) is the analytic function (of s) defined by 



r 00 °° / 



k 2 a 2 
At 



2 2 

" m 



?P(s; a; m) = a 

./o 

Making a change of variable i 1— » a 2 t, we have 

*P(s; a; m) = a 1+2s / t^ 1 V exp I -to 
^ fe=i v 

For s > 1/2, Jacobi inversion formula (|3.7p gives 

qj(s; a; m) =a 1+2s / i 8 " 1 £ e - taW i ~ + V^rt ^ 

~' fc 1= l I fc=-c 



dt. 



fe 2 \ 



4t / 



dt. 



dt. 



When a is small, we can use the Taylor expansion of e ta m to get 

-2s / y 



^P(s; a; m) 



-r(s) + v^to" 1 " 2 ^ s 



^ / *i \ j POO ^ 

,-=o J - Jo k=l 



-2 s 



-r(s) + Vrnn _1 ~ a T s 



2v^a 1+2 ^ 



j\ (2tt) 2 J+ 2s + 1 



T r + - C(2s + 2i + l), 



where £(s) is the Riemann zeta function. One can check that the formula on 
the right hand side of the last equality define an analytic function of s. Putting 
s = — (d + l)/2 and taking derivative with respect to a, we find that if d is odd, 



T(a;L 2 , ...,L d ;m) 



=2ir d ' 2 a- d - 1 



£ UlV^WV^r ( j + 1 - u - ) C (2j - d) 



Li=2 J j6NU{0} 
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and when d is even, 

T(a;L 2 , . . .,L d ;m) 



=27r d / 2 a- d - 1 



d+l ' 
2 , 



i=2 



,<i+l 



2 d +V 



j"GNU{0} 

^4 



(2^-W^r j + i-^ C(2j-d) 



(-I)t m d 



Here is t ne logarithm derivative of the gamma function, and for s < 0, T(s/2- 
1)C( S ) is understood as 



r 



s + 2 



7T 5 ' 2 / 1 — S 



-sT 



C(l- S ). 



When <i = 1,2, 3, using £(2) = 7r 2 /6 and £(4) = 7r 4 /90, we have respectively 

~i \ 7r 7rm 2 , _ n , 

S(a;m) ~ — loga + O{a ); 



and 



T(a;L 2 ,i3;m) 



30a 4 



rC(3) 



m 2 L 2 L 3 
24a 2 



m 2 £ 2 
47ra 



32 



+ O(a ); 



-L 2 i 3 log a + O(a ). 
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